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ABSTRACT 

The evolution and distribution of the angular momentum of dark matter (DM) halos have been dis¬ 
cussed in several studies over the past decades. In particular, the idea arose that angular momentum 
conservation should allow to infer the total angular momentum of the entire DM halo from measur¬ 
ing the angular momentum of the baryonic component, which is populating the center of the halo, 
especially for disk galaxies. To test this idea and to understand the connection between the angular 
momentum of the DM halo and its galaxy, we use a state-of-the-art, hydrodynamical cosmological 
simulation taken from the set of Magneticum Pathfinder simulations. Thanks to the inclusion of the 
relevant physical processes, the improved underlying numerical methods, and high spatial resolution, 
we successfully produce populations of spheroidal and disk galaxies self-consistently. Thus, we are 
able to study the dependence of galactic properties on their morphology. We find that (1) the specific 
angular momentum of stars in disk and spheroidal galaxies as a function of their stellar mass com¬ 
pares well with observational results; (2) the specific angular momentum of the stars in disk galaxies 
is slightly smaller compared to the specific angular momentum of the cold gas, in good agreement 
with observations; (3) simulations including the baryonic component show a dichotomy in the specific 
stellar angular momentum distribution when splitting the galaxies according to their morphological 
type (this dichotomy can also be seen in the spin parameter, where disk galaxies populate halos with 
slightly larger spin compared to spheroidal galaxies); (4) disk galaxies preferentially populate halos 
in which the angular momentum vector of the DM component in the central part shows a better 
alignment to the angular momentum vector of the entire halo; and (5) the specific angular momentum 
of the cold gas in disk galaxies is approximately 40% smaller than the specific angular momentum of 
the total DM halo and shows a significant scatter. 

Subject headings: dark matter - galaxies: evolution - galaxies: formation galaxies: halos - hydro¬ 
dynamics - methods: numerical 


1. INTRODUCTION 


The question of how the different types of galax¬ 
ies formed, and, in particular, which properties deter¬ 
mine a galaxy to become a disk or spheroidal, has been 
a matter of debate since the discovery that the Uni¬ 
verse is full of distant stellar islands of different mor¬ 
phology. For elliptical galaxies the following formation 
scenarios are discussed: from observations of ongoing 
mergers between galaxies of similar masses we know 
that spheroids can form in a major merger event be¬ 
tween (spiral) galaxies . This has been supported by 
many simulations (e . g.. lToomrelll977t !Whitdll978ln~979t 

iBarnes fe Hernaulstl 119961 : INaab et al.l 1200611 . However, 


this scenario is not sufficient to explain certain observa¬ 
tional properties, especially for the most massive galax¬ 
ies. Thus, an alternative formation scenario through a 
series of mu l tiple minor merger events was proposed by 
iMeza et all fl2003h and establ i shed bv. e.g.. INaab et al.1 
( 2009f). Gonzalez-Garcfa et all (f2009h . lOser et al.l ( 2010D . 
an d l.Tohansson et al.1 (|2012t h 
For the formation scenarios of spiral galaxies, the de¬ 
tails of their formation processes are less well known, but 


all the different channels discussed in the literature over 
the past decades are connected to the detailed buildup of 
the angular momentum and how the gaseous component 
transports intrinsic angular momentum into the central 
part of the halo, where the galaxy assembles. Initially, 
the dark matter (DM) halo and the infalling gas have 
identical angular momenta, but during the formation of 
the galaxy the gas cools and condenses in the center of 
the halo to form a disk. Assuming that the angular mo¬ 
mentum of the gas is conserved during this process, the 
angular momentum of the halo of a disk galaxy can be 
estimated directly f rom the angular momentum of the 
disk of the ga laxy (|Fall fe Efstathioiil 119801 iFalll 119831 
IMo et al.lH998lh 

In hierarchical scenarios of structure formation, struc¬ 
tures form through the gathering of c lumps owing to 
the gravitational force dPecblesI 119931 and references 
therein). The DM collapses at high redshifts into small 
objects, which grow into larger objects through merging 
and finally build halos. Those halos are not completely 
spherically symmetric owing to tidal torques induced by 
neighboring halos. The baryons condense in the cen- 
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ters of these DM halos and form the first protogalaxies 
(IWhite fc ReesllT978l) . Since the baryons and the DM 
originally gained the same a mount of angu l ar momentum 
throu gh these tidal torques ((Peebles 196f| iDoroshkevichl 
I1970T) and the angular momentum of the gas should be 
conserved during the collapse, the disk is expected to 
have a simila r specific angular mome ntum to that of its 
hosting halo (iFa.ll fe Efstathioulll980D . 

In contrast to observations, numerical Al-body simu¬ 
lations have the advantage that structures can be fol¬ 
lowed through time, enabling the detailed understand¬ 
ing of the early stages of galaxy formation. Until 
recently, simulations using traditional smoothed par¬ 
ticle hydrodynamics (SPH) codes suffered from the 
so-called angular momentum problem, where the ob¬ 
jects became too small compared to observations be¬ 
cause the gas had lost too much angular momen- 


turn Navarro & Benz 

19911 INavarro Kr. Whitdl1994 

iNavarro & Steinmetzl 1997 

). iSiiacki et al.l (]2012l) showed 


that SPH simulations overestimate the number of el¬ 
liptical galaxies owing to the lower amount of mixing, 
and therefore causing this spurious “angular momentum 
crisis” of the baryons. On the opposite side, adaptive 
mesh refinement (AMR) or moving mesh simulations 
tend to overpredict the amou nt of disk-like galaxies i n 
the absence of any feedback (|Scannanieco et all 120121) . 
As simulations become better in resolution and also in¬ 
clude feedback from stars, supernovae (SNe), and ac¬ 
tive galactic nuclei (AGNs), the gas is prevented from 
cooling too soon. Hence, early star formation is sup¬ 
pressed and the associated loss of angular momentum 
can be minimized. There are several studies that in¬ 
vestigate the influence of star formation and the associ¬ 
ated SN feedback on the formation of galactic disks (e.g.. 
Brook et al.ll2004 lOkamo to c t al.ll2005UGoycrnato et al'l 

2007t iScannanieco et al.ll2008HZavala et al.ll2008D . They 
find that strong feedback at early times leads to the for¬ 
mation of more realistic disk galaxies. A more recent 
study on the effect of stellar feedback on the angular mo¬ 
mentum is presented in llJbler et al.l ( 2014 1, where they 
found that strong feedback favors disk formation. A re¬ 
cent, det ailed summary o f disk galaxy simulations can be 
found in lMurante et al.l (120151) . 

One parameter to describe the rotation of a system 
is the so-c a lled sp in parameter A, which was introduced 
bv iPeeblc^ (1 9691 1971 1 and has since been investigated 
in several studies. With A it is possible to measure the 
degree of rotational support of the total halo. It has the 
advantage that it is only weakly de pending on the halo’s 
mass and its internal substructures dBarnes fc Efstathioul 
119871 ). The connection of this parameter with galaxy 
formation has been studied by many authors. 

At f irst, DM-o nly simulations were employed. 
iBullock et al.l (120011) introduced a modified version of the 
A-parameter by defining the energy of the halo via its cir¬ 
cular velocity. In addition, they studied the alignment 
of the angular momentum of the inner and outer halo 
parts. They found that most halos were well aligned 
but that abou t 10% of the halos showed a misalignment. 
lAubert et al.l (120041) investigated the alignment between 
the inner spin of halos and the angular momentum of 
the outer halo and also the alignment between the in¬ 
ner spin and the angular momentum of the inflowing 


material. The misalignment angle of the angular mo¬ 
mentum vectors at different r adii w as studied in more 
detail by iBailin fe Steinmetzl (120051) . who showed that 
with increasing separation of the radii the misa li gnment 
of the vectors increases. In another study, iMaccio et al.1 
( 20081 ) focused on the spin parameter’s dependency on 
the mass and the cos mology and found no correlations. 
iTrowland et al.1 (120131) studied the connection of the halo 
spin with the large-scale structure. They confirmed that 
the spin parameter is not mass dependent at low redshift 
but found a tendency to smaller spins at higher redshifts. 

Other authors included nonradiative gas in their anal¬ 
ysis. The misalignment of the gas and DM angular mo¬ 
mentum vectors and their spi n parameters at r edshif t 
z = 3 were investigated bv Ivan den Bosch et al.l (120021) . 
They found the angular momentum vectors to be mis¬ 
aligned by a median angle of about 30°. However, the 
overall distributions of the spin parameters of the gas and 
the DM were found to be very similar. An other detailed 
study bv lSharma fe Steinmetzl (120051 ) and iSharma et al.1 
(120121) found that the spin parameter of the gas com¬ 
ponent is on average higher than that of the DM, and 
they reported a misalignment of the angular momentum 
o f the gas and the D M of about 20°. 

I Chen et all (120031 ) compared the spin parameter and 
the misalignment angles between the DM and gas an¬ 
gular momentum vectors obtained in simulations, which 
include radiative cooling. This enables a splitting of the 
gas into a cold and a hot component. In their nonra¬ 
diative model they confirmed that the spin parameter of 
the gas component has higher spin than the DM, while in 
their simulations with cooling the two components had 
approximately the same spin. The misalignments of the 
global angular momentum vectors were 22.8° and 25° for 
the two differen t cooling mode l s and 23.5° for the non¬ 
radiative case. iStewart et al.1 ( 2011 .. 120131 ) focused on 
the specific angular momentum and the spin parameter 
in relation to the cold/hot mode accretion by following 
the evolution of individual halos over cosmic time. They 
found that the spin of cold gas was the highest compared 
to all other components and ab out 3—5 times high e r than 
that of the DM component. iScannapieco et all ;21)1)9; 
investigated the evolution of eight individual halos from 
zoom simulations and did not find a correlation between 
the spin parameter and the morphology of the galaxy. 
However, they found a correlation between the morphol¬ 
ogy and the specific angular momentum: for disks, the 
specific angular momentum can be higher than that of 
the host halo, while spheroids tend to have lower spe¬ 
cific angular momentum. In addition, they saw a mis- 
align ment betwe e n ste llar disks and infalling cold gas. 
Also. lSales et al.l ( 2012 ) reported no correlation between 
the spin parameter and the galaxy type and concluded 
that disks are predominantly formed in halos where the 
freshly accreted gas has similar angular momentum to 
that of earlier accretion, whereas spheroids tend to form 
in halos where gas streams in along misaligned cold flows. 
The angular mo mentu m prop e rties of the inflowing gas 
were studied by iPichon et al.l (j201lf ). They found that 
the angular momentum of the cold dense gas is well 
aligned with the angular momentum of the DM halo, 
in contrast to the hot diffuse gas. Additionally, an in¬ 
crease of the advected specific angular momentum of the 
gas component with cosmic time was noticed. They pro- 
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posed a scenario in which the angular momentum of the 
galaxies is fed by the collimated cold flows that are am¬ 
plified with time and make the disks larger. iKimm et al.l 
( 20Il f) studied the different behavior of the gas and DM 
specific angular momentum by following the evolution 
of a Milky-Way-type galaxy over cosmic time. They also 
found that the gas has higher speci fic ang u lar momen tum 
and spin parameter than the DM. lHahn et al.1 (120101) in¬ 
vestigated a sample of about 100 galactic disks and their 
alignment with their host halo at three different redshifts. 
Both the stellar and gas disks had a median misalignment 
angle of about 49° with respect to the hosting DM halo 
at z = 0. 

Some works focused explicitly on the effect that 
baryons have on the DM halo by running parallel DM- 
only simulations. Br^anmt^il. ( 120131) . who extracted ha¬ 
los from OWLS (jSchave et al.l l201(Tl . did not find a de¬ 
pendency of the spin parameter on mass, redshift, or cos¬ 
mology in their DM-only runs. In the simulations that in¬ 
cluded baryons and strong feedback, the overall spin was 
found to be affected only very little, while the baryons 
had a noticeable influence on the inner halo, independent 
of the feedback strength. iBett et all ( 2010 ) investigated 
the specific angular momentum and the misalignment of 
the galaxy with its halo. They obtained a median mis¬ 
alignment angle of about 25° for the DM-only runs and 
about 30° for the run including baryons. The baryons 
were found to spin up th e inne r region of the halo. 

Recently, I Welker et ah (120141 ) studied the alignment of 
the galaxy spins with their su rrou nding filament s, using 
the Horizon-AGN simulation dDubois et al.|[2014l) . They 
find that halos experiencing major mergers often lower 
the spin, while in general minor mergers can increase the 
amount of angular momentum. If a halo does not un¬ 
dergo any mergers but only smooth accretion, the spin 
of the galaxy increases with time, in contrast to that of 
its hosting DM halo. Since the gas streams and clumps 
in general move alon g the filaments, t he g alaxies realign 
with their filaments. iDanovich ct al.l (120151) have investi¬ 
gated the buildup of the angular momentum in galaxies, 
using a sample of 29 resimulated galaxies at redshifts 
from z = 4 to 1.5. Overall the spin of the cold gas was 
about three times higher t han that of t he DM halo, in 
line with previous studies. iGenel et al.l (120151) explored 
the importance of stellar and AGN feedback for the evo¬ 
lution of the angular momentum of the galaxies in the 
Illustris simulations, finding the angular momentum of 
their simulated galaxies to be in good agreement with 
observations for a chosen feedback similar to the one used 
in this paper. 

Observations indicate that the morphology of galax¬ 
ies is strongly influenced by t he relation between mass 
and angular momentum (see iFalll 11983 ). The angular 
momentum of disk galaxies was found to be about six 
ti mes hig her th an tha t of t he ellip ticals of equ al mass. 
iRomanowskv fe Falll (|2012l ) and iFall fc Romanowskvi 
(120131) revisited and extended this work, analyzing 67 
spiral and 40 early-type galaxies. They found that lentic¬ 
ular (SO) galaxies he between spiral and elliptical galax¬ 
ies in the so-called M sta r — jstar-plane. The bulges of 
spiral galaxies follow a simila r relation because they be¬ 
have like “mini-ellipticals.” lObreschkow fc Glazebrookl 
(120141) used data from high-precision measurements of 
16 nearby spiral galaxies to calculate the specific angu¬ 


lar momentum of the gas and the stars. They confirmed 
observationally that the mass and angular momentum 
are stron gl y correlated with the morp hology of galax¬ 
ies. iHcrnandez fc Cervantes-Sodil ( 2006 ) calculated the 
A-parameter for two galaxy samples with a total of 337 
observed spirals and found that with decreasing bulge- 
to-disk ratio the spirals have increa sing A-values. In a 
further studv lHernandez et all (120071 ) investigated a sam¬ 
ple of 11,597 spiral and elliptical galaxies from the Sloan 
Digital Sky Survey (SDSS) and found that ellipticals on 
average have lower A-values than spiral galaxies. 

Recently, a new generation of cosmological, hy- 
drod ynamical simu l ations , e.g., the MassiveBlack- 
II faha ndai et al.1 I2015D , Magneticum Pathfinder 
dHirschmann et all 120141). Illustris dVo gclsbcr ger et all 
120141). Horizon-AG N ( Dubois et al.ll2014D . and EAGLE 
(jSchave et al.ll2015l ) simulations, have been employed to 
follow the evolution of structures in the Universe. A new 
aspect of such simulations is that for the first time rea¬ 
sonable galaxy morphologies can be associated with the 
alaxies formed in t hose simul ations, e.g., Horizon -AGN 
Dubois ct ah . 2014 ). Illustris (jTorrev et al.l 12015 )). EA¬ 


GLE (iSchave cWd. 2015;), and Magneticum Pathfinder 
(iRernus et al.1 120*151 ). In this work we will analyze 
simulations from the set of Magneticum Pathfinder 
simulationtQ (K. Dolag et ah, in preparation), which 
are introduced in Section 2, and investigate how the 
baryonic component influences the morphology of the 
galaxy. In Section 3 we introduce the formulae and 
investigate the angular momentum, the spin parameter, 
and the alignments of halos. In Section 4 we show the 
kinematical split-up of galaxies, which allows us to clas¬ 
sify the galaxies in the simulation using the circularity 
parameter e in Section 5. We then examine the angular 
momentum of the different components of spheroids 
and disks. Furthermore, we study the alignment of the 
angular momentum vectors of the baryonic and the DM 
component. Additionally, we analyze the differences in 
the spin parameter A of spheroidal and disk galaxies. 


2. THE MAGNETICUM PATHFINDER SIMULATIONS 

In order to study the properties of galaxies in a sta¬ 
tistically relevant manner, we need both a large sample 
size and high enough resolution to resolve the morphol¬ 
ogy and underlying physics of galaxies. Since the newest 
generation of cosmological simulations can achieve both, 
they are a valuable tool for this study. We take the 
galaxies for our studies from the Magneticum Pathfinder 
simulations, which are a set of cosmological hydrody- 
namical simulations with different volumes and resolu¬ 
tions (K. Dolag et al., in preparation). The simula¬ 
tions were performed with an extended version of the 
Wbody/SPH code GA DGET-3, which is an updated 
versio n of GADGET-2 (iSprineel et al.l l2001bt iSprineell 
2005'). It includes various updates in the formulation of 
SPH regarding t he treatment of the viscosity and the 
used kernels (see iDolag et al.ll2005f iDonnert et al.1120131 : 
iBeck et aI1l2015t) . 

It also allows a treatment of radiative cooling, heating 
from a uniform time-dependent ultraviolet (UV) back¬ 
ground, and star formation with the associated feedback 
processes. The latter is based on a subresolution model 


1 www.magneticum.org 
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for the multiphase structur e of the interstellar medium 
(ISpringel fe Hcrnauistll2003r) . Radiative cooling rates are 
computed following the same procedure presented by 
IWiersma et al.1 ( 2009 1. We account for the presence of 
the cosmic microwave background (CMB) and of UV/X- 
ray backgrou nd radiation from q uasar s and galaxies, as 
computed by lHaardt fe Madaul (1200 IK . The contribu¬ 
tions to cooling from each one of 11 elements (H, He, 
C, N, O, Ne, Mg, Si, S, Ca, Fe) have been precomputed 
using the pub licly availab le CLOUDY photoionization 
code (jFerland et all 119981 1 for an optically thin gas in 
(photo) ionization equilibrium. 

In the multiphas e mo del for star formation 
dSpringel fe Hernauistl 120031 1. the ISM is treated as 
a two-phase medium where clouds of cold gas form from 
cooling of hot gas and are embedded in the hot gas phase 
assuming pressure equilibrium whenever gas particles 
are above a given threshold density. The hot gas within 
the multiphase model is heated by SNe and can evapo¬ 
rate the cold clouds. A certain fraction of massive stars 
(10%) is assumed to explode as SNe II. The released 
energy by SNe II (10 51 erg) is modeled to trigger galactic 
winds with a mass loading rate being proportional to 
the star formation rate (SFR) to obtain a resulting 
wind velocity of w w i n d = 350km/s. Our simulations also 
inc lude a detailed model of chemical evolution according 
to iTornatore et all ( 2007 1. Metals are produced by 
SNe II, by SNe la, and by intermediate- and low-mass 
stars in the asymptotic giant branch (AGB). Metals 
and energy are released by stars of different mass by 
properly accounting for mass-de pendent li fetimes (wi th 
a lifetime function according to iPadovani fc Matteuccil 
199, Ifl. the metal li city-d ependent stellar yields by 
Wooslev fe WeaveH (1199511 for SNe II, the yields by 
van den Hoek fe Groenewegenl (1 1997H for AGB stars, 
and the yields bv iThielemann et al.l (12003T) for SNe la. 
Stars of different mass are initially distributed acco rding 
to a Chabrier initial mass function (lGhabrieiIl2003 f. 

Most importantly, our simulations also include a pre¬ 
scription for black hole (BH) growth and for feed- 
back from AGNs b ased on the model p r esente d in 
iSpringel et al.1 (|2005f) and iDi Matteo et all (|2005h . in¬ 
cl uding the same modifications as in the study of 
iFabian et al.l (2010 :1 and some new, minor changes. 

As for star formation, the accretion onto BHs and the 
associated feedback adopt a subresolution model. BHs 
are represented by collisionless “sink particles” that can 
grow in mass by accreting gas from their environments, 
or by merging with other BHs. The gas accretion rate 
M. is est imated using the Bondi-Hoyle-Lyttleton approx¬ 
imation (iHovle fe LvttletorJ 1 1 9.39t iBondi fe Hovlelll944 
IBondilH952li~ 

„> _ 4tt G 2 M?f hoost p ^ 

( C 2 +Z, 2 ) 3 /^ ’ [i) 

where p and c s are the density and the sound speed of 
the surrounding (ISM) gas, respectively, /boost is a boost 
factor for the density, which typically is set to 100 and v 
is the velocity of the BH relative to the surrounding gas. 
The BH accretion is always limited to the Eddington 
rate (maximum possible accretion for balance between 
inward-directed gravitational force and outward-directed 
radiation pressure): M. = min(M., M e dd)- Note that 


the detailed accretion flows onto the BHs are unresolved, 
and thus we can only capture BH growth due to the 
larger-scale gas distribution, which is resolved. Once the 
accretion rate is computed for each BH particle, the mass 
continuously grows. To model the loss of this gas from 
the gas particles, a stochastic criterion is used to select 
the surrounding gas particles to be removed. Unlike in 
ISpringel et al.l (120051 1. in which a selected gas particle 
contributes with all its mass, we included the possibility 
for a gas particle to lose only a slice of its mass, which 
corresponds to 1/4 of its original mass. In this way, each 
gas particle can contribute with up to four ‘generations’ 
of BH accretion events, thus providing a more continuous 
description of the accretion process. 

The radiated luminosity L r is related to the BH ac¬ 
cretion rate by L r = e r M.c 2 , where e r is the radia¬ 
tive efficiency, for which we adopt a fixed value of 
0.1 (standardly assumed for a radiatively efficient ac- 
cretio n disk onto a nonrapi dlv spinning BH accord¬ 
ing tolSh akura fc Sunvaevlll973l ' see also ISoringell 120051 : 
IDi Matteo et al.l 120051 1. We assume that a fraction ef 
of the radiated energy is thermally coupled to the sur¬ 
rounding gas so that Ef = e r et M.c 2 is the rate of the 
energy feedback; ef is a f ree parameter and typic ally set 
to 0.1 (see discussion in iSteinborn et all ( 2015 11. The 
energy is distributed kernel weighted to the surround¬ 
ing gas particles in an SPH-like manner. Additionally, 
we incorpora t ed th e feedback prescription according to 
IFabian et al.l (I2010H : we account for a tran sition from a 
quasa r- to a radio-mode feedback (see also iSiiacki et al.l 
120011 ) whenever the accretion rate falls below an Ed¬ 
dington ratio of / e dd ; = M r /A(edd < 10 -2 . During the 
radio-mode feedback we assume a 4 times larger feed¬ 
back efficiency than in the quasar mode. This way, we 
want to account for massive BHs, which are radiatively 
inefficient (having low accretion rates), but which are ef¬ 
ficient in heating the ICM by inflating hot bubbles in 
correspondence to the ter mination of AGNjets. Note 
that we also, in contrast to ISpringel et al.l (j2005l) . mod¬ 
ify the mass growth of the BH by taking into account 
the feedback, e.g., AM. = (1 — rj r )M m At. Additionally, 
we introduced some technical modifications of the orig- 
in al implementation, for w hich readers can find details 
in iHirschmann et ahl (|2014T) , where we also demonstrate 
that the bulk properties of the AGN population within 
the simulation are similar to observed AGN properties. 

The simulati on addition ally follows thermal conduc¬ 
tion, similar to iDolag et al.l (l2004fl . but with a ch oice of 
1/20 of the classical Spitzer value (iSpitzerl Tl962h . The 
choice of a suppression value significantly below 1/3 can 
be justified by comparison with full MHD simulations in¬ 
cluding an anisot ropic treatment o f thermal conduction 
(see discussion in lArth et al.l ( 2014 /1. 

The initial conditions are using a standard ACDM 
cosmology with parameters according to the seven-year 
results of the Wilkinson Micro wave Anisotropy Probe 
(WMAP7] iKomatsu et aLlfeOllf l. The Hubble param¬ 
eter is h = 0.704, and the density parameters for matter, 
dark energy, and baryons are Dm = 0.272, Da = 0.728, 
and Df, = 0.0451, respectively. We use a normalization 
of the fluctuation amplitude at 8 Mpc of as = 0.809 and 
also include the effects of baryonic acoustic oscillations. 

The Magneticum Pathfinder simulations have al- 
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ready been successfully used in a wide range of nu¬ 
merical studies, showing good agreement with obser¬ 
vational results for the pressure profiles of the in- 
tracluster medium ( Planck Collaboration et all 1 201. 'll 


iMcDonald et, all 12014), for the prope rt ies of the AGN 
popu lation (|Hirschmann et all 12014 iSteinborn et all 
I2015D . and for th e dynamical properties of massive 
spheroidal galaxies (|Remus et al.ll2013l I2015D . 

In this work we mainly used a medium-sized (48 
Mpc//i) 3 cosmological box at the uhr resolution level, 
which initially contains a total of 2 • 576 3 particles (DM 
and gas) with masses of TOdm = 3.6 • 10 7 M & /h and 
m gas = 7.3 ■ lO 6 M 0 //i, having a gravitational softening 
length of 1.4 kpc /h for DM and gas particles and 0.7 
kpc /h for star particles. Additionally, we performed a 
DM-only reference run, where we kept exactly the same 
initial conditions, e.g., the original gas particles were 
treated as collisionless DM particles. 

To identify s ubhalos we used a version of SUBFIND 
(|Springel et al.l l2001af) . adapt ed to treat the baryonic 
component ( Dolag et al.l 120091) . SUBFIND detects ha¬ 
los based on a s tandard Friends-of-Friends algorithm 
dDavis et al.lll985 ) and self-bound subhalos around local 
density peaks within the main halos. The virial radius 
of halos is evaluated accordi ng to the density contrast 
based on the top-hat model ( Eke et all 119961 ) . In fur¬ 
ther post-processing steps we then extract the particle 
data for all halos and compute additional properties for 
the different components and within different radii, using 
the full (thermo) dynamical state of the different parti¬ 
cle species. In this study we mainly present results at 
redshifts z = 2, z = 1, z = 0.5, and z = 0.1. 


3. PROPERTIES OF HALOS 

From our data set we extract all halos with a virial 
mass above 5 • 10 u Mq. At a redshift of z = 2 there are 
396 halos, at z = 1 there are 606 halos, at z = 0.5 we find 
629 halos, and at 2 = 0.1 there are 622 halos. The lower 
limit is chosen to obtain a sample of halos that contain 
a significant amount of stellar mass and for which the 
resolution is sufficient to resolve the inner stellar, and 
gas structures. We transform the positions and velocities 
of all gas, stellar and DM particles into the frame of 
the halo, where the position given by SUBFIND is the 
position of the minimum of the potential and the velocity 
of the host halo is the mass-weighted mean velocity of all 
particles belonging to the halo. 

Exemplary, Figure Q] shows six galaxies within the cho¬ 
sen mass range demonstrating that the galaxies in the 
simulations look like observed disk and spheroidal galax¬ 
ies. The middle panel is a picture of the full cosmological 
box at redshift z = 0.5. 


3.1. The Angular Momentum 

The mass and the angular momentum of galax¬ 
ies are observed to be clo sely correlated wit h their 
morphology (e.g., iFalll 119831: iRomanowskv fc Fall 12012 : 

IQbreschkow fc Olazebrookl 120141) ' In a closed system 

without external forces, the angular momentum is a con¬ 
served quantity. However, in the context of galaxy for¬ 
mation and the interaction between collapsing objects 
and the large-scale structure, the assumption of the con¬ 
servation is not necessarily fulfilled for a single galaxy. 


The total angular momentum of a galactic halo is given 
by 

J = ^2 ( ^2 r M x V M ) ! ( 2 ) 

k \i£N k ) 


where k are the different particle types of our simula¬ 
tion (gas, stars, and DM) and N/. is their corresponding 
particle number with the loop index i. 

In our simulations, each particle carries its own mass. 
The initial mass is different for gas, star, and DM parti¬ 
cles. Later on, the mass of individual gas particles varies 
owing to mass losses during star formation. Thus, we re¬ 
move the mass dependence and use the specific angular 
momentum 


jfc = 


Yi m k,i r k,i X Vfc,» 

Yi m k,i 


(3) 


where k are the species of matter, as above. Therefore, 
we firstly calculate the angular momentum of each par¬ 
ticle of a species. Afterward, we sum over all individual 
particles and divide by the total mass of the correspond¬ 
ing species to obtain the absolute value. 


3.2. The Spin Parameter A 

In the following section we want to study the dimen - 
sionless A-parameter. As defined bvlPeeblesI (119691 . 11971D 
and adopted by, e.g., I Mo et all (119981) . the general A- 
parameter, which is used for the total halo, is given by 


GM 5 / 2 ’ 


(4) 


where E = —GM 2 /2 R V1I is the total energy of the halo. 
This total spin parameter can only be used considering 
all matter inside the halo. To evaluate the different com¬ 
ponents (gas, stars, and DM), th is param eter needs to be 
modified. We follow [Bullock et all (12001D . who defined a 
component-wise spin parameter as follows: 


A'(r) 


J( r ) 

V / 2M(r)U c irc(i’)r ’ 


(5) 


where J(r) is the angular momentum, M(r) the mass, 
and Wire = y/GM(r)/r the circular velocity within a ra¬ 
dius r. Another advantage of A '(r) is that it can be used 
for the calculation at different radii. When calculated 
over the entire viria0 radius, A'(r) ~ A. For simplic¬ 
ity we will drop the prime for the remaining part of our 
study. For the evaluation of the different components A 
can be expressed in terms of t he specif i c angu lar momen¬ 
tum, as done bv lvan den Bosch et al.l (120021) : 


Afc 


Jk 


\/2R v irWir 


( 6 ) 


where k stands for the different components. 

Independent of the definition, the distribution of the 
A-parameter can be fitted by a lognormal distribution of 
the following form: 


P(A) 


1 

Xy/2n a 



ln 2 ( A/A 0 )\ 

2 < t 2 y 


(7) 


2 Formally, if the true energy E is used, this is equal for a trun¬ 
cated, isothermal sphere. Ho wever, for a NFW h alo the correc¬ 
tion term is of order unity, see lBullock et all (1200II ) and references 
therein 
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Fig. 1.— Main panel: complete cosmological box at redshift z = 0.5. The upper panels show exemplary spiral galaxies, and the lower 
panels show spheroidal galaxies. We show two random projection directions for each of the galaxies. In all panels we color-code the stellar 
particles by their cosmological formation epoch, where the upper color bars represent the age of the stars from old to young and the lower 
color bars represent the gas temperature from cold to hot. We also marked the position of the halos within the cosmological box. 
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Fig. 2.— Left panels: histograms for the A-parameter for the 
different components within it v i r ; the DM (black) is the dominant 
component, the stars (yellow) peak at lower values, whereas the gas 
(turquoise) is distributed around higher values. While Ao stays rel¬ 
atively constant for the DM and stellar component, it increases for 
the gas with decreasing redshift. Right panels: the gas component 
splits into hot (red) and cold (blue) gas. Both components spin 
up with cosmic time, whereas for the cold gas this happens faster. 
Overplotted are the best-fit lognormal distributions (smooth lines), 
for which the parameters are reported in Table [l] 

with the fit parameter Ao, which is about the median 
value, and the standard deviation a. 

The left panels of Figure [2] show the histograms of 
the Afc-distribution in linear bins of all halos in the se¬ 
lected mass range for stars (yellow), DM (black), and gas 
(turquoise) from redshift z = 2 (upper panels) to z = 0.1 
(lower panels). The histograms and fit curves (smooth 
lines) are each normalized to the number of halos. The 
distribution of spin parameters of the stellar component 
has significantly lower values compared to the distribu¬ 
tion of the spin para meter for DM, in agreem ent with the 
results presented bv iDanovich et all 120151 1 . This is due 
to the fact that within the halo the stars are more con¬ 
centrated toward the center, whereas the spin of the DM 
component is dominated by the outer part of the halo, 
where most of the angular momentum of the DM com¬ 
ponent resides, as also shown in the right panel of Figure 
1221 in Appendix . In addition, major mergers result in 
a reduction of the spec ific angular momentum, as shown 
in IWelker et al.l (120141) . On the contrary, the gas com- 


TABLE 1 

The Ao Fit Value for the Different Components at 
Different Redshifts 


Redshift 

-^stars 

-^DM 

-^gas 

^hot 

-^cold 

2 

0.023 

0.043 

0.061 

0.053 

0.074 

1 

0.026 

0.046 

0.078 

0.069 

0.097 

0.5 

0.021 

0.042 

0.085 

0.078 

0.107 

0.1 

0.020 

0.042 

0.090 

0.083 

0.123 


ponent always has a spin parameter distribution shifted 
towar d larger values, in agre e ment w ith more rece n t stud - 


ies bvISharma & Steinmetz 

(20051). Kinnn et al.1 

201 if), 

Sharma et al.l (20121). and 

Danovich et al. (12015 

), but 


-—- v , - i - | - I I - w — 

in con trast to previous s tudies by Ivan den Bosch et all 
(120021 ) and IChen et ahl (|2003f ) . who found nearly the 


same spin distributions for gas and DM. The larger spin 
values for the gas reflect the continuous transport of the 
larger angular momentum from the outer parts into the 
center due to gas cooling. While this leads even to a spin¬ 
ning up of the gas component with time, the distributions 
for the stars and DM remain relatively constant. This 
spin-up of the gas component might be caused by the 
continuously accreted cold gas, which brings in large an¬ 
gular momentum fr om farther away along the filaments 
dPichon et al.ll201ll) . When dividing the gas into hot and 
cold phases, we note that the hot gas (red) has lower val¬ 
ues than the cold gas (blue), which has a long tail to high 
A-values. This is in contrast to IChen et al.l (2003), who 
found that the spin parameter for the hot gas is higher 
than that of the cold gas. On the other h and, our results 
agree well with more recent studies by IDanovich et al.1 
(|2015l) . Th is trend for the cold gas to high values was 
also seen bv IStewart et all 120111 ) , who have calculated A- 
values for the cold gas and found values around 0.1 — 0.2. 
For a better overview we have listed the fit value of the 
A-distribution for the different components in Tabic [I] 

3.3. Alignments 

In this subsection we want to study the orientation of 
the angular momentum vectors of the different compo¬ 
nents in comparison to each other within the innermost 
region of the halo. For the calculation of the innermost 
region of the halo we only take particles within the inner 
10% of the virial radius (R v i r ) into account, which corre¬ 
sponds roughly to the size of the galaxies. We chose this 
radius as an approximate medium value, since smaller 
disks have a size of about 5% while extended gaseous 
disks can reach out up to 20% of the virial radius. 

The angle between two vectors of the species i and j 
is calculated by 


We then bin the angles from 0 to 180 degrees in 18 
bins of a size of 10° each and count the number of halos 
within each bin. For the plots the number is normalized 
to unity. 

As shown in Figure [3l the distribution of the an¬ 
gles between the angular momentum vectors of gas and 
stars (green dotted) generally shows a good alignment at 
higher redshifts and becomes less aligned with decreas¬ 
ing redshift. This is in contrast to the alignment be¬ 
tween the DM and stellar components (yellow dashed), 
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Fig. 3. — Angle between the total angular momentum vectors 
of the different components within the innermost 10% of fJ v i r as 
indicated in the plots for different redshifts. The gray dot-dashed 
line is expected for a random distribution of the angles. At higher 
redshift the gas and stellar components (green dotted) are well 
aligned, and their alignment gets worse with decreasing redshift. 
In contrast, the alignment of the DM and the stars (yellow dashed) 
gets better with decreasing redshift. 

which shows a better alignment with decreasing red¬ 
shift. The angles between the DM and the gas com¬ 
ponent (turquoise solid) are an intermediate case, and 
their distribution does not change with redshift. For 
comparison, a random distribution of alignment angles 
is shown as a gray dot-dashed line in Figure [3l demon¬ 
strating that none of the angle distributions found in the 
simulations are of random nature. The difference be¬ 
tween the behavior and especially the evolution of the 
alignment of the gas and the stellar component can be 
explained by the fact that the angular momentum of the 
gas reflects the freshly accreted material (similar to the 
stars at high redshift), while the stars at low redshift 
reflect the overall formation history, similar to the DM 
content. These trends are in qualitative agreement with 
results obtain ed from si mulations at very high redshift 
(z > 9) by iBiffi fc Maiol (120131) , who also found the an¬ 
gular momentum of the gas component in the center to 
reflect the recent accretion history. 

4. KINEMATICAL SPLIT-UP 

In this section we investigate properties of our galaxies 
as a function of their stellar mass and angular momen¬ 
tum. We therefore calculate the angular momentum J of 
the stars within the innermost 10% of the virial radius. 
Furthermore, we ignore all particles within the inner¬ 
most 1% of the virial radius because of their potentially 
unknown contributions caused by bulge components or 
numerical resolution effects. We rotate the positions and 
velocities of all particles such that their ^-components are 
aligned with the angular momentum vector J of the stel¬ 
lar component. In the same manner we produce another 
data set that is rotated such that the ^-components are 
aligned with the angular momentum vector J of the gas 
component. These data are used for the e gas -distribution 
because the gas circularity calculated with the data ro¬ 


Fig. 4.— Stellar mass in the inner 10% of the virial radius vs. 
the specific angular momentum of the stars for all halos at redshift 
z = 0.1. The mean circularity £ s tar of the halos is col o r-cod ed. We 
adopted the scaling relations of IRomanowskv fc Falll (12012 b (their 
Figure 2), where the blue line stands for pure disks, the turquoise 
one for SO, and the yellow one for pure bulges. 

tated according to the angular momentum of the stars 
can fail to detect a gas disk that is misaligned with the 
stellar disk and vice versa. 

4.1. The Circularity Parameter e 

Since in the Magneticum simulations disk and 
spheroidal galaxies are formed, we need to categorize 
them for our analysis according to their morphology. In 
order to distinguish between different types of galaxies, 
we use the circularity p arameter e. T he e -parameter 
was first introduced by lAbadi ct al.l (120031 ) as Ej = 
J z /JrWciE). For our s tudy we use the definition of 
iScannapieco et all (120081) . which is given by 


Jcirc r Vcirc 

where j z is the ^-component of the specific angular mo¬ 
mentum of an individual particle and j c ire is the expected 
specific angular momentum of this particle assuming a 
circular orbit with radius r around the hal o center o f 
mass, with an orbital velocity of V r c i rc (r) = \/GM(r)/r. 

We compute the circularity e for every individual par¬ 
ticle between 1% and 10% of the virial radius. This is 
done for the stellar and the gas component, where we 
use the data that were rotated according to the angu¬ 
lar momentum of the stars or the gas, respectively. To 
obtain the circularity distribution for our selected halos, 
we compute the fractions /(e) of particles within equal- 
distant bins, using a bin size of Ae = 0.1. In a dispersion- 
dominated system there is usually a broad peak in the 
distribution at e ~ 0, while in a rotation-supported sys¬ 
tem there is usually a broad peak at £ ~ 1. 

To test the hypothesis that the mass and the angu¬ 
lar momentum of galaxies are the most important in¬ 
gredients in their formation h istory and the resulting 
morphology (c.g., iFalll 119831 IRomanowskv fc Falll l2012t 
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Fig. 5. — Averaged ^-distributions of the stellar component for 
four redshifts. Thereby, each e bin is averaged over all halos that 
lie in the corresponding 6-value bin. We see a clear transition of 
the galaxy types from rotational-supported systems to dispersion- 
dominated systems at all redshifts. 

lObrcschkow fc Olazebrooll I2014H . we plot our galaxies 
at z = 0.1 in the stellar mass vs. stellar angular mo¬ 
mentum plane in Figure 21 Thereto we take all star 
particles within the innermost 10% into account. We 
color-code the simulated galaxies according to the abso¬ 
lute 0 value of the mean of their stellar circularity pa¬ 
rameter. A dopting the c lassifi cation diagram of galaxies 
by iRomanowskv fc Falll (j2012f) (their Figure 2), we over¬ 
plot lines where the specific angular momentum follows 
the relation j oc M a with a ss 2/3. We plot different 
lines with distance A6 = 0.25. As can clearly be seen in 
Figured! the mean circularity parameter is following this 
relation. Inspired by this result, we classify our galaxies 
according to what we will in the following refer to as the 
“6-value”: 

ft = >^(»)-5 log '“(vf)- (10) 

which is the y-intercept of the linear relation f(x) = 
ax + b in the log-log of the stellar-mass-specific-angular- 
momentum plane._ 

According to IRomanowskv fc FalU (I2012D . galaxies 
with 6-values close to —4 (blue line) are expected to be 
disks, while 6 « — 5 indicates pure spheroidals (yellow). 
In our simulations we even find galaxies with 6-values 
down to 6 = —6.25 for the galaxies with the smallest 
specific stellar angular momentum. 

Figure [5] shows the averaged e-distribution of the stel¬ 
lar component at four redshifts, colored according to the 
different 6-value bins shown in Figure |U Each e bin is 
averaged over the halos that reside in the chosen 6-value 
bin. There is a clear transition between the galaxies with 
different dynamics, from rotational-supported ((e) = 1) 

3 We take the absolute value, since there is a tiny fraction of 
halos with small 6-values for which e can have a slightly negative 
value owing to the different weightings in the definitions of the 
specific angular momentum (eq. @) and the circularity (eq. ©)■ 



Fig. 6. — Top: mean circularity e of the star component av¬ 
eraged over all halos in the corresponding 6-value bin. There is 
a clear trend for decreasing mean e when moving from the pure 
disks (high 6-value) down to the pure bulges (smaller 6-value) for 
all four redshifts. Middle: same as the top panel, but for the e of 
the gas component. Bottom: The mass fraction of cold gas with 
respect to the stars, both within the inner 10% of the virial radius, 
averaged over all halos in the corresponding bin. At high redshift 
(pink solid line) there is much more cold gas, even in halos with a 
smaller 6-value. 

to dispersion-supported ((e) = 0) systems, reflected by 
their 6-value from the M star -jstar-plane: While galaxies 
with 6 ps —4 clearly have most stars around e = 1, and 
thus are dominated by rotation, galaxies with 6 sa — 6 
peak around e = 0, indicating that there is no significant 
rotation. This is true for all redshifts. However, we see 
a slight trend with redshift at intermediate 6-values: for 
example, the distribution of the halos in a 6-value bin of 
—4.75 (green squares) at z = 2 has the higher of the two 
peaks around e = 1, whereas at 2 = 0.1 it only peaks 
around e = 0. It is also interesting that the galaxies in 
the 6-value bin of —4.5 (turquoise stars) at redshift z = 2 
have a dominant rotationally ordered component, which 
becomes more dispersion dominated with decreasing red¬ 
shift. We also note that the distributions around £ = 0 
become slightly more distinct at lower redshifts. This in¬ 
dicates that there is an evolution of the different galaxy 
types with cosmic time along the M s t a r-jstar-pla-ne. 

This additional evolutionary trend gets more evident 
in Figure [HI which shows £ s t ar (top), £ gas (middle), and 
the mass fraction of the cold gas with respect to the 
stellar mass (bottom) versus the 6-value, with spheroidal 
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systems at the left and disk-like systems on the right. 
Thereby, each of these properties is averaged over the ha¬ 
los in the corresponding 6-value bin. We can clearly see 
that the average £ sta r increases with increasing 6-value at 
all four redshifts. Additionally, the average stellar circu¬ 
larity is generally larger at higher redshifts for all types 
of galaxies, albeit this trend is stronger for galaxies with 
larger 6-values. This clearly shows that disk galaxies at 
higher redshifts had less prominent central bulges that 
would shift the value toward (e) = 0. For (e gas ) we do 
not find any clear trend. The bottom panel of Figure [6] 
shows a similar trend for the mean cold gas fraction. The 
spheroidal systems at present day with small 6-values 
have only low amounts of gas, compared to the stellar 
mass, usually below 10%, whereas disks have gas frac¬ 
tions of 20% or more. For higher redshifts this fraction 
increases successively, and the strongest evolution is vis¬ 
ible between redshifts z = 2 (pink solid line) and 2 = 1 . 
Below 2 = 1 there is only a mild but continuoutQ evo¬ 
lution, but at 2 = 2 even spheroidals have gas fractions 
of 20% or more, while disk galaxies can contain more 
than 40% gas compared to their stellar content. For the 
largest 6-values the galaxies can even be dominated by 
gas, i.e., the gas fraction is larger than 50%. Thus, we 
conclude also that galaxies that are spheroidal systems 
have a nonnegligible cold gas fraction at high redshifts. 

4.2. The X-parameter 

We want to understand whether this classification of 
the galaxies according to their position on the M star -j s t a r- 
plane has an effect on the spin parameter A. In Figure 
[3 we show the distribution of the total A-parameter of 
the whole halo content (see Equation |4|), exemplary for 
redshift 2 = 0.1. We divide the halos according to their 
6-value. The first interesting thing to note is that there 
is a double peak in the A-distribution for galaxies with 
disk-like kinematics (upper left panel). This could indi¬ 
cate that there are different formation channels for disk¬ 
li ke systems, e.g., gas-ri ch major m ergers, as discussed 
i nlSnringel & Hcrnauistl (|2005D and lSchlachtberger. D.P.I 
(120141) . The peak at higher A-values becomes smaller as 
we move to the right and seems to disappear at a 6-value 
of —4.75. Galaxies with 6-values below —5.25 are hosted 
in halos with significantly smaller A-values compared to 
the overall distribution (black dashed line). Therefore, 
on average there seems to be a connection between the 
galaxy types and the spin parameter A of the hosting 
halos, reflected in a continuous transition of the spin pa¬ 
rameter distribution of the host halo with the 6-value 
obtained from the M s t ar -j s tar-plane. 

4.3. Alignments of the Central Components 

To verify the dynamical connection between the differ¬ 
ent baryonic components within the galaxies, we show 
in Figure [8] the distributions of the alignment between 
the central angular momentum of the stellar and the gas 
component. As before, we classify our galaxies according 
to the 6-value. By moving from high (upper left) to low 

4 At z = 1 (magenta dashed line) the value for the galaxies with 
the highest 6-value exceeds the scale, since there are six galaxies 
that on average have a cold gas fraction of 2.53 with respect to the 
stars. 


(lower right) 6-values we immediately see that the align¬ 
ment between the angles is good for the disk galaxies 
(large 6-values) and gets increasingly more random for 
decreasing 6-values. This is indicated by the gray dotted 
line, which again illustrates a random distribution (as 
in Figure |3j) . We find a continuous transition from the 
disk-like to the bulge-dominated galaxies. It illustrates 
that in spheroidal systems the stars are significantly mis¬ 
aligned with an eventually present gaseous disk. This re- 
sult is su p ported by the findings of the ATLAS 3D -project 
(jCappellari et al.il201ll) that a kinematical misalignment 
of the g as componen t with respect to the stars is not 
unusual dDavis et al.ll201lD . 

5. CLASSIFICATION OF SIMULATED GALAXIES 

So far we have seen that there is a continuous tran¬ 
sition between the different types of galaxies within the 
Af s t ar -Jst a r-plane, orthogonal to j oc M a with a ~ 2/3. 
We now want to study the kinematical properties of our 
simulated galaxies depending on their classification as 
disk and spheroidal galaxies. 

5.1. Selection Criteria 

In our simulations, many spheroidal galaxies show ex¬ 
tended, ring-like structures of col d gas, in good agree¬ 
ment with recent observations (jSalim fe Richl 120101) . 
Therefore, including the circularity of the gas within 10% 
of the virial radius as a tracer of morphology can lead 
to a misinterpretation. Additionally, there can be huge 
uncertainties for galaxies that have almost no gas left. 
Using a criterion only based on the 6-value of the galax¬ 
ies also is not straightforward because of the existence 
of fast rotators among the spheroidal galaxies. Observa- 
tionally there seems to exist some overlap between the 
different galaxy types within the M s t ar -jstar-plane. In 
principal, the luminosity of the stars could be taken into 
account, since old stars that make up the bulge and the 
halo stars are not as luminous as young stars that build 
up the disk. Hence, in observations a galaxy with an old 
massive stellar bulge and an extended disk of young stars 
and gas is very likely to be classified as a spiral galaxy. 

However, here we will stick to a classification of galax¬ 
ies based on the circularity distribution (Equation [2]) of 
their stars, which allows us to capture rotationally sup¬ 
ported stellar disk structures or dispersion-dominated 
spheroidal structures. We combine this criterion with the 
mass fraction of the cold gas with respect to the stars, 
following our result from Figure [6] As before, circular¬ 
ity distributions are evaluated within 10% of the virial 
radius, excluding the central 1%, while we include the 
central 1% in our calculation of the cold gas fraction as 
the resolution is not important in this case. We use the 
results from the previous section to justify the thresh¬ 
old values used to select dispersion-dominated, gas-poor 
spheroidal and rotational-supported, gas-rich disk galax¬ 
ies to reflect counterparts of classical, observed elliptical, 
and spiral galaxies. 

As shown in Figure 0 there is a clear, bimodal be¬ 
havior of the epsilon distribution within different 6-value 
bins. We show in Appendix in detail the cumulative dis¬ 
tributions that allow us to define the proper thresholds 
bracketing the transition regions. These thresholds are 
then applied to the circularity distributions of the indi¬ 
vidual galaxies, which of course in general show a more 
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Fig. 7.— Distribution of the total A-parameter for the halos according to their 6-values at z = 0.1. There is a transition from the disks 
(upper left) to the bulges (lower right), i.e., from rotation-dominated systems (higher A-values) to dispersion-dominated systems (lower 
A-values). Interestingly, the halos with the highest b -values show a dichotomy in their distribution. 



Fig. 8.— Angle between the angular momentum vector of the gas and that of the stars within the innermost 10% of the virial radius 
according to their 6-values at redshift z = 0.1. There is a clear transition between the galaxy types. The higher the b- value, the better 
aligned the stellar and gas components are. When going to lower b- values, i.e., moving down on the M s tar-istar-pl&ne to the bulges, the 
angles of the two components seem to be randomly distributed (compare gray dotted sinusoidal distribution). 
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complex shape than the averaged distributions. The left 
panel of Figure [20] in Appendix shows the cold gas mass 
fraction of the galaxies classified only by their circularity 
Estarj illustrating our choice of including a gas criterion. 
In short, we are using the following selection criteria: 

• We classify a galaxy as a spheroidal galaxy if the 
majority of particles is in a close interval around 
the origin (i.e., /(—0.4 < e < 0.4) > 0.6). This 
percentage cut is adapted to the redshift. For the 
details see the left panels of Figure [T9] in Appendix 

In addition, it has to have a cold gas mass fraction 
with respect to the stars lower than 0.35 at z = 2, 
0.2 at z = 1, 0.125 at 2 = 0.5, and 0.065 at z = 0.1. 
This criterion is chosen such that the fraction val¬ 
ues satisfy the linear function f(z) = 0.15-2 + 0.05. 
The upper limit for the spheroidal ga laxies at red- 
shift z = 0.1 seems plausible, since lYoung et all 
(12014ft found that massive elliptical galaxies (of the 
red sequence) have mass fractions of HI and H 2 , 
compared to the stars, up to 6% and 1%, respec¬ 
tively. 

• We classify a galaxy as a disk galaxy if the majority 
of particles are off-centered from the origin (i.e., 
/(0.7 < £ < 3) > 0.4; for the dependence on the 
redshift see right panels of Figure [TO] in Appendix 
)• 

Additionally, the constraint on the mass fraction of 
the cold gas is such that it has to be higher than 
0.5 at z = 2, 0.35 at z = 1, 0.275 at z = 0.5, and 
0.215 at z = 0.1. This criterion is chosen such that 
the values of the mass fraction satisfy the linear 
function f(z) = 0.15 • z + 0.2. 

• All remaining galaxies, which fulfill neither of the 
two criteria, are classified as “others.” 

The total number of galaxies above the mass cut of 
5 • lO n M 0 selected from our simulation at redshift 2 = 
0.1 is 622 (for other redshifts see Table [2]). According 
to the above criteria, 64 galaxies are classified as clas¬ 
sical disks and 110 as classical spheroids, while 448 are 
classified as “others.” Such unclassified galaxies include 
merging objects, bulge-dominated spirals, spheroids with 
extended gas disks, barred galaxies, or irregular objects. 

TABLE 2 

Number of Halos in the Selected Mass Range 


Redshift 

N 

-^spheroid 

-^Vdisk 

2 

396 

34 

89 

1 

606 

87 

73 

0.5 

629 

146 

59 

0.1 

622 

110 

64 


Note. -— Total number of halos N in the selected mass range, 
the number of spheroids -/V sp heroid and disks TVdisk (classified by 
the circularity of the stellar component and the mass fraction of 
the cold gas with respect to the stars) at different redshifts. 

In order to identify parts of the unclassified objects, the 
classification could be refined in such a way that we ad¬ 
ditionally divide into spheroids that have a gaseous star¬ 
forming disk or galaxies that have a large stellar disk with 


only little gas (SO galaxies). At lower redshifts our clas¬ 
sification becomes increasingly difficult since there are a 
lot of galaxies that have a very dominant stellar bulge 
but, on the other hand, also possess an extended gaseous 
disk containing many young stars. 

Although the criterion for the gas mass fraction is rel¬ 
atively arbitrary and may overestimate the number of 
spheroids, we checked that changing this selection crite¬ 
rion does not change the results qualitatively, although 
adding/removing galaxies from the spheroidal sample. 
At higher redshifts an additional difficulty might be that 
the transition of galaxy types is less strict, as seen in Fig- 
ure[5] At these redshifts, remaining dynamical signatures 
of the current formation process will be more pronounced 
and individual distributions of the circularity parameter 
might be more complex. Thus, a clear assignment might 
fail. 

5.2. Comparison of the Simulated Stellar Specific 
Angular Momentum with Observations 

After having extracted a set of spheroidal and disk 
galaxies, we can evaluate the relation between the spe¬ 
cific angular momentum and the mass of the stellar com¬ 
ponent and compare them with observations, as shown 
in Figure [9[ To compare w ith observations presented 
by (Fall Sz Romanowskvl (12013ft (hereafter FR13) we cal¬ 
culate the mass and specific angular momentum for all 
stellar particles within 10 R \/2 instead of 10%i? v ir, where 
R 1/2 is the radius that contains half the stellar mass of 
the galaxy and roughly corresponds to observed effective 
radii. This is done to account for the fact that, for a 
given mass, dis k galaxies have a larger effective radius 
than spheroids (jShen et al.l [2003ft . and to better resem¬ 
ble the radius ranges studied in FR13. The left panel 
shows the M s tar-jstar-plane for our galaxies classified as 
spheroids (red circles), including the 23 elliptical galaxies 
(orange circles) presented in FR13 (their Figure 2). In 
the right panel the same is shown for our disk galaxies 
(blue diamonds) in comparison to the 57 spiral galaxies 
(purple diamonds) from FR13. In general, our simulated 
galaxies are in good agreement with the observations. 
This result has already been shown for a subset of our 
galaxies at z = 0 with a rat her crude classification cri¬ 
terion in I Remus et all (12015ft . and we find an even bet¬ 
ter agreement with observations with our more adv anced 
classification scheme. Recently, IGenel et al.l ([2015ft have 
shown a comparison of the galaxies in the Illustris sim¬ 
ulation with those observations, which are also in good 
agreement when a similar feedback mechanism is used. 

To understand the impact of the choice of radius used 
to calculate the angular momentum and the stellar mass, 
we included the values for all galaxies evaluated within 
10% of the virial radius, as previously shown in Figured 
For the spheroidal galaxies the effect is much stronger, 
and we clearly see that the angular momentum is larger 
for larger radii. The radius dependence is shown in more 
detail in the right panel of Figure [22] in Appendix , where 
we show that the angular momentum of spheroids in¬ 
creases with radius and that there is a significant con¬ 
tribution to the specific angular momentum from large 
radii. For disk galaxies the exact radius is less important, 
as most of the angular momentum is contained within the 
disk at relatively small radii. The same holds true for the 
5-value, which more strongly depends on the considered 
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Fig. 9. — Galaxies within the stellar-mass-specific-angular-momentum plane at 2 : = 0.1. Gray symbols represent all galaxies extracted 
from the simulations (as in Figure|4|, where angular momentum and stellar mass are measured within 10% of the virial radius. The left panel 
shows the simulated galaxies classified as spheroids (red circles) compared to data from observed elliptical galaxies bv IFall &; Romanowskvl 
(120131) (orange circles, denoted by “FR13”). The right panel shows disks in our simulation (blue diamonds) compared to observational data 
of spiral galaxies by FR13 (purple diamonds). For the disks and the spheroids, we evaluated the mass and the specific angular momentum 
within 10 R\/ 2 , which more closely resembles the radii studied by FR13. The colored lines correspond to the b -values (see Figure [4]l. 


radius for spheroids than for disks, which can be seen in 
Figure [23] in Appendix . 

Note that we compare here the specific angular mo¬ 
mentum directly measured from the total stellar com¬ 
ponent from the simulations (Equation [3j), while the 
observations are inferred from the projected measure¬ 
ments. We also do not resolve galaxies with stellar 
masses smaller than ss 10 1 o Mq, while the observations 
include some objects with smaller masses. 

5.3. The Gas and Stellar Specific Angular Momentum 
in Simulations and Observations 

We also compare the specific angular momentum of 
the gas with that of the stars for all classified galaxies in 
our simulations. Figure [TO] shows the relations for all four 
redshifts as indicated in the plots. At redshift z = 0.1 we 
include the observational data taken from th e THINGS 
(The HI Nearby Galaxy Survey, IWalter et al.ll2008h sam¬ 
ple, which consists of 16 spiral galaxies in the local Uni¬ 
verse, for whi ch su rfac e dens ities of stars and cold gas are 
available (|Lerov et all 1200 811. Specific angular momenta 
for the gas a nd stellar components of tho s e g alaxies were 
presented bv IQbreschkow fe GlazebrooU (12014H . 

For the disk galaxies, we find that the specific angular 
momentum of all stars in the central region is slightly 
smaller than that of the gas (blue diamonds), in agree¬ 
ment with the observations (purple symbols). This, most 
likely, originates from the fact that the specific angular 
momentum of the gas is constantly replenished by freshly 
accreted gas, which transports larger angular momentum 
from the outer parts of the halo to the center. This be¬ 
comes more evident by looking at the newly formed stars 
(turquoise diamonds) in our simulations. We consider a 
star to be young if its formation happened not more than 
Az = 0.05 x 1/(1 + z) ago at a given redshift z, which 
corresponds to a stellar age of 0.2Gyr at z = 2, 0.4Gyr 
at z = 1, 0.5Gyr at z = 0.5, and 0.7Gyr at z = 0.1. If 
we only take the young stars into account, we find al¬ 
most an equality (dotted line) with the specific angular 
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Fig. 10. — Specific angular momentum of the gas against the 
specific angular momentum of stars, both within 10% of the virial 
radius for galaxies which are classified as disks (blue diamonds) at 
four redshifts as indicated in the plots. Additionally, we show the 
relation also considering only young stars (turquoise diamonds). 
At z = 0.1 we over plot observa tional data points calculated by 
IQbreschkow &; Glaze brook! (120141) (purple pentagons). The values 
agree well with the observations, and we note an overall spin-up 
with cosmic time. 

momentum of the gas. At a lower redshift (lower right 
panel) the values of the specific angular momentum of 
the gas are slightly larger than at higher redshift (up¬ 
per left panel). This behavior reflects the spin-up of the 
cold gas with time, as already seen before. Hence, also 
the young stars have higher specific angular momentum. 
This appears in a slight trend for a separation of the 
young stars and the total stellar component, which have 
lower specific angular momentum than the gas. 
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The spheroids (red circles in Figure ITOl) . however, have 
a significantly lower specific stellar angular momentum 
compared to that of their gas. This suggests that, es¬ 
pecially at high z, most gas in spheroids originates from 
the accretion of material from large radii (cold streams or 
infalling substructures) and hence transports the higher 
angular momentum from the outer parts into the center. 

We can clearly see that the spheroids and the disk 
galaxies show a different behavior in the relation between 
the angular momentum of the gas and stars. We also 
have seen that the gas gains angular momentum over 
time, especially in disk galaxies. In disk galaxies, stars 
have slightly smaller specific angular momentum than 
the gas, and they also show a mild difference between the 
specific angular momentum of stellar components with 
different ages. Here the young stars have slightly larger 
specific angular momentum, basically reflecting the an¬ 
gular momentum of the gas from which they form. 

All in all, the results from the simulations fit well with 
the observ ational dat a. In particular, we are in agree¬ 
ment with iFalll (119831) , who find that the specific angu¬ 
lar momentum of the galaxies increases with the disk-to- 
bulge ratio for a given mass. 

We note that the evaluation of this quantity for 
spheroids especially at lower redshifts might be fault- 
prone since, owing to the selection criterion, there is only 
a small amount of cold gas within the chosen radial range. 

5.4. Comparison of the Gas and DM Specific Angular 
Momenta 

Finally, we compare the specific angular momenta of 
DM and gas. In particular, we are interested in the 
scaling relations of the corresponding individual specific 
angular momenta of baryonic and nonbaryonic matter. 
Previous studies (e.g., lFalllll983l : IMo et al.Hl998fl suggest 
an equality between the angular momentum of the total 
DM halo and that of the central gas component of disk 
galaxies. 

The relation between the specific angular momentum 
of the cold gas of the galaxy, which resides within the 
innermost 10% of the virial radius, and that of the total 
DM halo is shown in Figure QT] for different redshifts. In 
general, we find the specific angular momentum of the 
DM to exceed (by a factor of « 2) the specific angular 
momentum of the cold gas. Interestingly, the cold gas 
in disk and spheroidal galaxies does behave in the same 
way. This is due to the fact that cold gas always settles 
in disk-like structures, even in elli ptical galaxies (e.g., 
iLees et all 119911 : lYoung et all 1201 111 . We plot all disks 
above the redshift-dependent M co id/M sta r-cut (blue dia¬ 
monds). Additionally, we fit a line parallel to the 1:1 re¬ 
lation (blue dotted) to the data, in order to see the offset. 
For all shown redshifts we find values of the y-intercept 
between 0.33 and 0.41 in log space, which is between 
2.14 kpc km/s and 2.57 kpc km/s in normal space. How¬ 
ever, the scatter around this fit curve is very large. The 
yellow diamonds show the value for one sim ul ated galaxy 
at different redshifts as presented in iKimm et all ( 20 Ilf ) . 
Their disk galaxy has only slightly higher specific angular 
momentum in the gas component compared to the DM 
halo, but sits within the distribution of values we find 
in our simulation for different galaxies. Our values are 
also b roadly in line with the results of iDanovich et ahl 
(l2015f) . who find that the spin of the disk is comparable 
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Fig. 11. — Specific angular momentum of the cold gas within 
10 % of the virial radius against the specific angular momentum of 
DM within the entire virial radius. The blue diamonds show the 
disks and the red circles the spheroidal galaxies, where the size of 
the circles reflects the fraction of the cold gas mass with respect 
to the stellar mass. The blue dotted line shows the fit for the disk 
galaxies. The 1:1 relation is represented by th e black dashe d line. 
For comparison, we add the results obtained bv IKimm et al.l (1201 il l 
(yellow diamonds). 


to that of the DM halo. The size of the symbols of the 
spheroids (red circles) is plotted according to their frac¬ 
tion of cold gas mass with respect to the stellar mass. 
For completeness we plot all spheroids with a cold gas 
mass fraction lower than 35% for all redshifts. At higher 
redshift most spheroids have a large amount of cold gas. 
At the lowest redshift there are still some objects that 
were classified as spheroids with the £ s tar-criterion but 
have a high cold gas mass fraction. Interestingly, there 
are some spheroids with high specific cold gas angular 
momentum, which might be due to individual infalling 
clumps and small substructures, as these are spheroidal 
galaxies that have extremely small cold gas fractions. 
However, in such cases we cannot relate the specific an¬ 
gular momentum of these individual structures to the 
specific angular momentum of the halo. 


5.5. Misalignment Angles 

As seen before, another indicator for the different evo¬ 
lutionary states reflected by the stars, DM, and gas 
are the misalignment angles between these three com¬ 
ponents. We now investigate their behavior by focusing 
on our selected disk and spheroidal galaxies. 


5.5.1. The Angle between Gas and Stars 

At first, we investigate whether we find a correlation 
between the morphological type and the angles between 
the angular momenta of the gas and the stars, both 
within the innermost 10% of the virial radius. Fig¬ 
ure m shows the angle for four different redshifts. The 
disk galaxies (blue dot-dashed histograms) have very well 
aligned gas and stellar angular momentum vectors with 
median values of 7.8° at redshifts z — 2 and z = 0.1 and 
median values of 6.3° and 7.6° at redshifts z = 1 and 
z = 0.5, respectively. This is in good agreement with 
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Fig. 12.— Angle between the angular momentum vector of gas 
and stars, both within the innermost 10% of the virial radius. At 
all redshifts the disks (blue dot-dashed) are well aligned. In con¬ 
trast, the spheroids (red dashed) are randomly distributed. The 
overall distribution (green solid) gets worse aligned with decreasing 
redshift. 


TABLE 3 

Median Misalignment Angles between the Baryons in the 
Center 


Redshift 

All Halos 

Disks 

Spheroids 

2 

13.4 

7.8 

37.4 

1 

24.4 

6.3 

56.5 

0.5 

35.8 

7.6 

65.5 

0.1 

37.6 

7.8 

60.4 


Note. — The median misalignment angles between the angular 
momentum vectors of gas and stars, both within the inner 10% of 
the virial radius (6> gas io-starsio)- 

lHahn etall ( 2010 ). who found median angles of about 
8° for their disks at z = 1 and a median of 7° at red- 
shifts ^ = 0.5 and z = 0. It demonstrates that the stellar 
and gaseous disks in disk galaxies are very well aligned, 
which is what we expect, since the stars form out of the 
gas and thus maintain the same orientation. 

The spheroids have a random distribution with a me¬ 
dian value of 60.4°, i.e., the gas and stellar components 
in spheroidal galaxies are often misaligned. In summary, 
the gas and star components of spheroids become less 
aligned with decreasing redshift, while there is no change 
for disk galaxies (see Table 0). The overall trend for all 
halos in the simulation shows the same behavior as the 
spheroids, in agreement with Figure [Sj 

5.5.2. The Angle between DM and Baryons 

Since we see a clear difference in the alignment between 
the angular momenta of the stars and gas for spheroids 
and disks, we want to test whether this is reflected in the 
relation between the angular momentum of the baryonic 
components and DM. The left panel of Figure [13] shows 
the misalignment angle between the angular momentum 
of the total DM halo and the gas within the inner 10% 
of the virial radius at redshift z = 0.1. Again, the gray 
dashed line is the expected distribution if the angles were 


spread randomly. We clearly see that the median mis¬ 
alignment found for all halos (green solid lines) is ps 49°, 
for the disks ps 45°, and for the sp heroids ~ 47°. This i s 
slightly larger than the results of iSharma et all 1 2012 ). 
who found a median misalignment angle of about 30° 
for the gas within the innermost 10% of the virial radius 
compared to the total ang ular momen t um v ector of all 
halos. On the other hand, lHahn et al.l (120101) compared 
the angular momentum of the gas component of the disk 
to the total angular momentum and obtained about 49° 
at z = 0, well in line with our results. 

We now want to see under which circumstances the 
orientation of the angular momentum vectors of the to¬ 
tal DM halo is reflected in that of the stellar component, 
since we have seen in Figure [12] that the angular momen¬ 
tum vectors of the gas and that of the stars are very well 
aligned in the inner part of the halo for disk galaxies and 
poorly aligned for spheroids. The angle between the an¬ 
gular momentum vectors of the DM and that of the stars 
within the virial radius is shown in the middle panel of 
Figure [T3l In general, there is an alignment of the two 
components. The disk galaxies are slightly better aligned 
with a median angle of ps 18° compared to the spheroids 
that have a median angle of ps 31°. As a median value for 
all halos we find « 24° at redshift z = 0.1 and a similar 
value for z = 1, namely, ps 22°, which is not shown here. 

The right panel of Figure [13] shows the angle between 
the angular momentum vector of the stars within the 
inner 10% of the virial radius and the DM within the 
entire virial radius. The distribution of the alignment 
angles looks similar to that of the gas (left panel of Figure 
m for all galaxies. The alignment seems poor, with 
median angles of ps 55° for all halos, ps 46° for disks, 
and ps 57° for sph eroidal galaxies. This agrees well with 
lHahn et al.l 1 2010 ), who calculated a median angle of the 
stellar component and the total halo conte nt of ab out 
ps 49° at z = 0 for disk galaxies, similar to iCroft et al.l 
(|2009f) reporting a media n angle of ps 44° at redshift 
2=1. iBett et ah! (120101) found slightly smaller angles 
of ps 34° for the alignment of the galaxies with respect 
to the hosting t otal DM ha l o. Th e galaxies at z = 1.2 
investigated bv iCodis et al.l (l2015f ). using the Horizon- 
AGN simulation, are slightly less aligned with their DM 
halos than our simulated g alaxies (see also the r ight panel 
of Figure[20]of Appendix ). iDeason et al.l (120111 1 reported 
that 41% of their disk galaxies had misalignment angles 
larger than 45° with their DM halo. In addition, they 
found that the disk galaxies are better aligned with the 
DM in the innermost 10% of the r^oo- Though we do 
not show this, here we obtained a median angle of ps 9° 
for the angular momentum vectors of the stellar and DM 
components in the innermost 10% for our disk galaxies 
at 2 = 0.1, and thus we can confir m this trend. It also 
agrees well with lHahn et al.l (12010 ). who find a median 
value of ps 15° for the angle between the DM and stellar 
components of the galactic disks. For the DM and the 
gaseous components of their disks they obtain a median 
value of ps 18°, which is in line with the median value of 
ps 12° for our disk galaxies. 

6 . THE HOST HALOS OF DIFFERENT GALAXY TYPES 

So far we have seen the intrinsic difference of the bary¬ 
onic components and how they reflect global halo prop¬ 
erties for galaxies classified as either disk or spheroidal 
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Fig. 13.— Left: angle between the total angular momentum vector of DM in R v ; r and the gas within the innermost 10% of it v ; r . The 
overall alignment is poor. Middle: The angle between the total angular momentum vectors of the DM and stars within R v \ r . The disk 
galaxies (blue dot-dashed lines) seem significantly better aligned than the spheroids (red dot-dashed lines). Right: angle between the total 
angular momentum vector of the DM in _R v ir and the stars within the innermost 10% of the virial radius. For all shown distributions the 
alignment is very poor, almost random. 
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Fig. 14.— Angle between the total angular momentum vector 
of the DM in i? v ; r and that of the DM within the innermost 10% 
of the virial radius. At redshifts z = 2, z = 0.5, and z = 0.1 we 
overplot the distribution of the DM-only run (black dotted). 


galaxies. We now want to study whether there are un¬ 
derlying differences in the DM halo contributing to the 
formation of these two classes of galaxies. 

6.1. The Alignment of the Host Halos with Their 
Centers 

One interesting question is if there are signatures of 
differences in the internal structure of the angular mo¬ 
mentum in the DM component for disk galaxies and 
spheroidal galaxies. Therefore, we compare the angular 
momentum of the DM within the whole halo with the 
one within 10% of the virial radius. Figure [14] shows the 
cumulative distributions of the angle between the DM 
angular momentum vector of the whole halo and that of 
the inner part of the halo. At z = 2, z = 0.5, and z = 0.1 
we include the corresponding distribution of the DM-only 


TABLE 4 

Median Misalignment Angles of the DM at Different Radii 


Redshift 

All Halos 

Disks 

Spheroids 

DMO 

2 

56.7 

54.8 

57.2 

66.6 

1 

54.9 

51.2 

61.3 

— 

0.5 

48.3 

53.1 

48.8 

53.5 

0.1 

47.4 

42.2 

50.6 

51.3 


Note. — The median misalignment angles between the angular 
momentum vectors of the DM within 10% and that of the total 
DM halo ($dm—dmio) f° r the hydrodynamical run as well as for 
the run with only DM (DMO). 

run 0- The alignment is generally better for disks than 
fo r spheroids (see also Table [4]), as already speculated 
in I Bullock et al.1 (120011) . Interestingly, the relative dif¬ 
ference between the alignments for disk and spheroidal 
galaxies is largest at z = 1, which corresponds to a typ¬ 
ical formation redshift of spiral galaxies. Overall, the 
misalignment grows with redshift for all halos. At red- 
shift z = 0.1 we calculate me dian values for all halos of 
~ 47°, which agrees well with lHahn et al.l ([20100 . report¬ 
ing a value of 45°. 

In contrast, in the DM-only run the values are slightly 
higher at all redshifts, with median misalignment an¬ 
gle of_~_51° at_z = 0.1. This tendency was also seen 
bv lBett et al.l (120101) . who found that in their run with 
baryons the vectors were slightly more aligned than in 
the DM-only case. The misalignment angle in our analy¬ 
sis is higher than their median angles of 15° —30° for their 
run with baryons. This could be due to the fact that we 
only consider the inner 10% instead of 25%, as i n their 
study. This is expected, since iBailin fc Steinmctd (|2005l ) 
found that the alignment becomes worse when the radii 
are further separated. 

A possible interpretation of the above findings could 
be that disk galaxies preferentially reside in halos, where 
the core is b etter ali gned w ith the outer parts of the 
halo (see also lBullock et akll200lD . It might well be that 
in such halos the angular momentum can be transported 
more effectively by the cooling of gas from the outer parts 
into the central parts. Another possibility could be that 

5 For technical reasons there were no data available at z = 1 
from the DM-only run. 
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disk galaxies survive merger events longer, when consec¬ 
utive infall is aligned with the angular momentum of the 
galaxy. Interestingly, spheroidal galaxies (especially at 
z ~ 1) show exactly the opposite behavior. The inner 
and outer parts of their DM halos are less aligned, in¬ 
dicating that major merging events are contributing to 
the buildup of spheroidal galaxies. This is in line with 
previous studies showing that major mergers with mis¬ 
aligned spins c an be responsible fo r angular momentum 
misalignments (l.Shajrna_et__aL|[20 Ill). 

In addition. IWelker et al.l if2014T ) proposed that 
anisotropic cold strea ms rea li gn th e galax y with its host¬ 
ing filament. However. iSales et~ahl (12012H suggested that 
disks form out of gas having similar angular momentum 
directions, which would favor the spherical hot accre¬ 
tion mode, while the accretion along cold flows that are 
mainly misaligned tends to build a spheroid. To answer 
this question in detail, further investigation is needed to 
trace back disk galaxies and see whether the primordial 
alignment causes the inflowing matter to become a spi¬ 
ral or whether it is the other way around, i.e., if the 
galaxy type over cosmic time induces the alignment. We 
suspect that the environment has a significant impact 
on the formation of disk galaxies. In less dense environ¬ 
ments, where the halos can evolve relatively undisturbed, 
the angular momentum of the galaxy preferably remains 
aligned with that of its host halo, and thus a disk can 
form. 

6.2. The Spin Parameter A of the Host Halos 

In the following section we finally return to the spin 
parameter A, evaluated for the whole halo. In Figure 
[l5l we show the distribution of the A-parameter for the 
total matter distribution within our halos for different 
redshifts. The histograms and lognormal fit curves (see 
Equation 0 are each normalized to the number of the 
halos. The distribution of A for all 396 halos at z = 2 
(upper left), 606 halos at z = 1 (upper right), 629 ha¬ 
los at z = 0.5 (lower left), and 621 halos at z = 0.1 
(lower right) is shown in green. We also plot the dis¬ 
tribution of the spheroidal (red dashed) and disk (blue 
dot-dashed) galaxies. The spheroids, which have only 
little rotation in the stellar component, tend to have 
lower A-values. The disk galaxies have their median at 
higher values. This bimodality of the two ga laxy types is 
seen a t all four redshifts, as already shown in lTcklu et al.l 
• 2013 :. The spheroids have always lower median values 
than the disk galaxies (see a lso Table[5]for the fittin g pa¬ 
rameters) . Observationally, iHernandez et al.l (I2007T) also 
found that the spiral galaxies on average have larger A- 
values than the ellipticals of their sample Howev e r, this 
was not seen in previous stu dies by ISales et id! (|2012T ) 
and iScannanieco et all (12009( 1 , who did not find a corre¬ 
lation between galaxy type and the spin parameter. 

The A-distribution for all halos seems to be relatively 
constant with time, i.e., independ e nt of redshift, which 
is in agreement with lPeirani et afl (l2004f) . 

To verify the statistical significance of the differences, 
we show in Figure [16] the cumulative distributions for 
the A-values, again splitting the halos according to the 
different galaxy types they host. This illustrates the dis¬ 
tances between the two distributions, where the curve 
for the halos of the spheroidal galaxies is always left of 
the curve for all halos and the curve of the disk galax- 



Fig. 15. — The A-distribution calculated with formula 0 for the 
different redshifts as indicated in the plots. Green takes all halos 
into account. The distributions of the spheroids are the red dashed 
histograms, and those of the disks are blue dot-dashed histograms. 
The curves are the fits given by equation q. At all redshifts a 
split-up of the two different galaxy types is present. The spheroids 
have their median at the lower A-values, while the disks have higher 
values. 


TABLE 5 

Calculated Median Values and the Fit Values of the 
A-Distributions 


Redshift 

All Halos 



2 

Amed = 

0.043, A 0 

= 0.043, 

a = 0.662 

1 

^med 

0.046, A 0 

= 0.047, 

a = 0.635 

0.5 

^med 

0.041, A 0 

= 0.042, 

er = 0.640 

0.1 

^med 

0.042, A 0 

= 0.043, 

a = 0.630 

Redshift 

Disks 




2 

^med 

0.058, A 0 

= 0.059, 

a = 0.543 

1 

^med 

0.060, Ao 

= 0.062, 

o- = 0.644 

0.5 

^med 

0.051, A 0 

= 0.050, 

a = 0.617 

0.1 

^med 

0.069, A 0 

= 0.064, 

a = 0.570 

Redshift 

Spheroids 



2 

^med 

0.028, A 0 

= 0.029, 

a = 0.497 

1 

^med 

0.042, A 0 

= 0.042, 

a = 0.633 

0.5 

^med 

0.037, A 0 

= 0.039, 

a = 0.620 

0.1 

^med 

0.034, A 0 

= 0.034, 

a = 0.546 


Note. — The calculated median values A me d and the fit values 
Ao and a of the A-distributions at different redshifts for our sample 
of galaxies. 

ies stays always to the right. The differences are more 
pronounced for small spins. To quantify the statistical 
significance of differences in the distributions for the disk 
and spheroidal galaxies, we apply a Kolmogorov-Smirnov 
(K-S) test. This test is applied to the two unbinned dis¬ 
tributions for the disks and the spheroids. Table [S] shows 
the calculated values for the maximum distance D and 
the probability. With the exception of the result^ at 
z = 0.5, this confirms that the spin distributions of the 
halos hosting disk galaxies are statistically significantly 

6 Although unlikely, the value obtained at z = 0.5 does not allow 
to exclude the same origin of the two distributions. 
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TABLE 6 

K-S Test for the A-Distributions 


Redshift 

D 

Probability 

2 

0.563 

1.41 -10- 7 

1 

0.319 

4.45 -lO” 4 

0.5 

0.237 

1.45 TO -2 

0.1 

0.421 

6.31 TO -7 


Note. — The maximum distance D and the significance level 
(probability) resulting from the K-S test that the A-distributions 
for disks and spheroids at different redshifts for our sample of galax¬ 
ies are from the same distribution. 



Fig. 16. — The A-parameter plotted against the normalized cu¬ 
mulative number of halos. Green takes all halos into account, red 
the spheroidal galaxies, and blue the disk galaxies. The distri¬ 
butions of the disks and the spheroids do not overlap (besides at 
z = 0.1 at the higher end). The disks are on the right side, i.e., at 
higher values, of the distribution for all halos, while the spheroids 
are on the left, i.e., at lower values. 


different from the halos hosting spheroidal galaxies. 

In this section we have seen that there is a statistical 
correlation between the morphological type and the over¬ 
all distribution of the spin parameter A. We thus con¬ 
clude that the total DM halo somehow “knows” about 
the morphology of the galaxy at its center. At all four 
considered redshifts the distributions for the spheroids 
have lower median A-values than those of the disks. Ac¬ 
cording to the K-S test, there is a strong indication that 
they do not originate from the same distribution. On the 
other hand, we have also found that there are spheroids 
with even higher A-values than disk galaxies. 

To verify that these differences are intrinsic to the halo 
and not caused by the contributions of the baryonic com¬ 
ponent to the spin parameter, we calculated the spin pa¬ 
rameter Adm for the DM component of the halos. Figure 
[IT] shows the AoM-distribution for the baryon run for the 
four redshifts as indicated in the plots. The green his¬ 
togram shows Adm for all halos, the red one stands for the 
spheroids, and the blue one for disk galaxies. Here again, 
most prominently visible at redshifts z = 2 and z = 0.1, 
there is a splitting of the two different galaxy types, the 
distribution of the spheroids peaks at lower values, and 



0 0.04 0.08 0.12 0.16 0.20 0.04 0.08 0.12 0.16 0.2 


Adm A dm 

Fig. 17.— The A-parameter for the DM component within R v i r ; 
also in the DM component of all halos (green) there is a broad 
distribution, which splits into spheroids (red) at lower values, while 
disks (blue) peak at higher values. 


TABLE 7 

K-S Test for the Adm-D^tributions 


Redshift 

D 

Probability 

^DMO 

Probability DMO 

2 

0.537 

5.92 10" 7 

0.590 

8.14 TO -s 

1 

0.319 

4.45 -10 -4 

— 

— 

0.5 

0.234 

1.65 TO" 2 

0.285 

2.3 TO -3 

0.1 

0.390 

4.60 TO- 6 

0.370 

9.84 TO” 5 


Note. — The maximum distance D and the significance 
level (probability) resulting from the K-S test that the Adm - 
distributions for disks and spheroids at different redshifts originate 
from the same distribution for the baryon run and the DM-only 
(DMO) run. 

that of the disks at higher values. We also perform a K-S 
test on the two distributions (see Table 0. The values 
are similar to that of the general A-distribution (Table 

0 . 

In the Appendix [C] we show that this split-up of the 
galaxy types is also reflected in the spin parameter A of 
the stellar component. 

In order to see whether the differences originate from 
differences within the DM halo or if they are caused by 
the interplay between the DM and the baryonic compo¬ 
nent, we finally compare the run with the baryons to the 
DM-only control run. We thereto cross-identify the ha¬ 
los of the DM-only run with those of the hydrodynamical 
run. We search for the corresponding halos in the DM 
run, allowing the center of the halo to be in a range of 200 
kpc around its position in the simulation with baryons, 
and additionally restrict to halo pairs that only differ by 
up to 30% in their virial masses. At redshift z = 2 we 
found 364 halos in the baryon run that match with the 
DM-only run, at z = 0.5 we could assign 609 halos, and 
at z = 0.1 we cross-matched 575 halos. We find that 
the AnM-values of each halo identified in both runs are 
very similar (see Figure [TT] of Appendix ). Since baryons 
have an effect on the DM, in the two runs the evolution 
of some of the halos can be quite different, as also dis- 
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cussed in lBett et al.l (120101) . Therefore, we cannot match 
every halo. 

Figure [TSl shows the A-distribution for the DM compo¬ 
nent at three different redshifts. The overall distribution 
is shown in black. The fit values Ao for the DM-only 
run are slightly lower than for the run with baryons. 
When we split the halos in the DM-only run according 
to their galaxy type assigned in the run with baryons, 
we find a similar split in the distributions as in the run 
with baryons. At redshift z = 2 we could use 76 disks 
and 33 spheroids, at z = 0.5 there are 55 disks and 143 
spheroids, and at z = 0.1 we cross-matched 53 disks and 
99 spheroids. The results of the K-S test (see Table 0 
show that these two distributions are unlikely to origi¬ 
nate from the same one. It is striking that even in the 
run without baryons, the split-up of the spin parame¬ 
ter for the different galaxy types is clearly visible. This 
suggests that the hosting DM halo and, connected with 
that, the formation history and environment play an im¬ 
portant role for the morphology of the resulting galaxy. 

7. DISCUSSION AND CONCLUSIONS 

We extracted between 400 and 630 halos with total 
halo masses above 5 • 10 11 Mg at four different redshifts 
from the hydrodynamical, cosmological state-of-the-art 
simulation Magneticum Pathfinder, which includes de¬ 
tailed treatment of star formation, chemical enrichment, 
and evolution of supermassive BHs. We investigated the 
distribution of the spin parameter A& of the different 
components (DM, stars, hot and cold gas) and the align¬ 
ment of the angular momentum vectors of those compo¬ 
nents within the entire halo, as well as within the central 
10% of the virial radius. To classify the galaxies accord¬ 
ing to their morphology, we rotate them such that their 
angular momentum vector is oriented along the 2 -axis. 
For this orientation we calculate the circularity parame¬ 
ter e, which allows a classification based on their circu¬ 
larity distribution. This allows us to define the subset 
of galaxies clearly identified as spheroidal or disk galax¬ 
ies and to compare their properties to observations. We 
additionally performed and analyzed a DM control sim¬ 
ulation to test the effect of the baryonic processes and 
the formation history on the angular momentum within 
the halos. We summarize our findings as follows: 

(i) For all our halos, the stellar component generally 
has a lower spin than the DM component, while the gas 
shows a significantly higher spin parameter, especially 
the cold gas component that dominates the overall spin 
parameter of the gas. While the distribution of the spin 
parameters of stars and DM does not show a signifi¬ 
cant evolution with time, the spin distribution of the gas 
component significantly evolves toward larger spin values 
with decreasing redshift. 

(ii) In general, the angular momentum vectors of the 
baryonic components and the DM are well aligned. There 
is an evolution within these alignments, where at high 
redshifts the alignment between gas and stars is better 
than their alignment with the DM component, while at 
low redshift stars and DM tend to be better aligned. 

(iii) When classifying galaxies according to their posi¬ 
tion within the stellar-mass-specific-angular-momentum 
plane, we demonstrate that various galaxy proper¬ 
ties show smooth transitions, as expected when go¬ 
ing from rotation-dominated systems to dispersion- 


dominated systems. This is most prominently found 
for the circularity distribution of stars and gas, the cold 
gas fraction, but also for the spin parameter distribution 
and the alignments of the angular momenta. Rotation- 
dominated galaxies have generally larger spin parame¬ 
ters, and the angular momentum vectors of the stars and 
gas components are aligned. On the contrary, dispersion- 
dominated galaxies generally have smaller spin parame¬ 
ters and show almost no alignment between the angular 
momentum vectors of stars and the (however small) gas 
components. 

(iv) Alternatively, when classifying our galaxies ac¬ 
cording to their circularity distribution in combina¬ 
tion with the cold gas content into disk galaxies and 
spheroidal galaxies, they populate clearly distinguish¬ 
able regions within the stellar-mass-specific-angular- 
momentum plane. Using that classification, we lose a 
significant number of galaxies that cannot be classified 
within this scheme. Nevertheless, for those galaxies that 
can be classified the results are in excellent agreement 
with observations for both spheroidal and disk galaxies. 

(v) In disk galaxies the specific angular momentum 
of the gas is slightly higher than that of the stars, 
which js in good agreement with recent observations 
(lObreschkow fe GlazebrooU2014D . This is due to the fact 
that the gas component has contributions from freshly 
accreted gas with higher specific angular momentum, 
which will be turned into stars later. This is also re¬ 
flected in our result that the specific angular momentum 
of the young stars is slightly larger compared to that of 
all stars. 

(vi) In general, the specific angular momentum of the 
total halo is higher than the specific angular momentum 
of the gas in the galaxy. For disk galaxies it accounts for 
roughly 43% of the value found for the total halo, with 
no significant redshift evolution. 

(vii) Overall, the angular momentum of the total halo 
is only weakly aligned with the angular momentum of 
the central part. Here the simulation including baryons 
shows slightly more alignment than the DM-only simula¬ 
tion. However, we found that in general the halos hosting 
disk galaxies (in comparison with halos with spheroidal 
galaxies at their centers) show a better alignment of the 
angular momentum vectors of the total halo and the 
central part for the DM component. This is most pro¬ 
nounced at redshift z ~ 1, where most disk galaxies are 
forming. 

(viii) The splitting of the galaxies into disk and 
spheroidal galaxies reveals also a dichotomy in general 
halo properties: the halos hosting disk galaxies have a 
slightly larger spin than the halos hosting galaxies clas¬ 
sified as spheroidal galaxies. This dichotomy is even 
reflected in the distribution of the spin parameters in 
the DM control run, where we cross-identified the halos 
which in the hydrodynamical simulation host galaxies of 
different types. This indicates that the formation history 
of the DM halo plays an important role for defining the 
morphology of the galaxies. 

Our results are based on the classification of galax¬ 
ies obtained from the circularity of the stellar compo¬ 
nent (e s tar) in combination with the fraction of cold gas 
with respect to the stellar mass within the central part 
of the halo, where the galaxies form. This allows us to 
select classical disk and spheroidal galaxies, which show 
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Fig. 18.— The A-parameter for the DM component for the DM-only run (black) at three redshifts. The red curves show the distribution 
for halos that were identified in the baryon run and classified as spheroids, while the blue curves show the halos classified as disks in the 
baryon run. There is a split-up of the galaxy types, which suggests that the morphology could be a result of the formation history. 


very distinct dynamical properties, reflected in their cir¬ 
cularity distribution, cold gas fractions, specific angu¬ 
lar momentum, and spin and angular momentum vector 
alignments. They therefore shed light on the main forma¬ 
tion mechanism of these galaxies within the cosmological 
framework. 

While our current classification focuses on a selection of 
galaxies whose properties resemble those of classical spi¬ 
ral and elliptical galaxies, the lion’s share of the galaxies 
in our simulation cannot be classified as poster child disk 
or spheroidal galaxies. This again reflects observational 
facts, since at present day many galaxies exhibit signs of 
distortions, ongoing merger events, peculiar structures, 
or other irregularities. For our unclassified galaxies, we 
clearly see that their properties show a smooth transi¬ 
tion between the poster child disks and spheroids, and 
to understand those transition processes, more complex 
classification schemes are needed. We suggest that those 
schemes should consider a differentiation between young 
and old stellar components, as well as the general distri¬ 
bution of the diverse gas phases in different parts of the 
galaxies. 

Additionally, when comparing to observations, more 
emphasis needs to be placed on how observational quan¬ 
tities are obtained and how this is mimicked when analyz¬ 
ing simulations. Classifications of galaxy morphologies in 


cosmological simulations also allow the interpretation of 
cosmological simulations in various new ways, as already 
proven in this work, not only for studying the forma¬ 
tion and evolution of galaxies but also for studying the 
relation of AGNs with their host galaxies, as well as for 
cosmological studies, where the bias of the measurements 
depends on the morphologies of the galaxies. Therefore, 
the good agreement of the intrinsic dynamical properties 
for our classified galaxies with respect to current obser¬ 
vations can be seen as a first step to promote cosmologi¬ 
cal, hyclrodynamical simulations for future cosmological 
studies. 
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APPENDIX 

APPENDIX A: DETAILS ON THE CLASSIFICATION 

In order to illustrate the choice of the cuts in £ sta r used for the classification of our galaxies, we show the cumulative 
circularity distribution for spheroids (left panels) and the ‘anti-’cumulative distribution for disk galaxies (right panels) 
in Figure 1191 For the spheroids we sum up the fractions from —£ to e of the distributions shown in Figure [5l 
f(x) = Yl-x fraction^). The cut is drawn between the halos of the 5-value-bin that show a clear behavior of 
spheroids at the corresponding redshift. The same is done for the disk galaxies, besides that here we sum up the 
fraction from e to 3, f(x) = f raction (e). 

In the left panels of Figurej20l we show that at high redshift there are many spheroidal galaxies (red histograms, 
classified with £ star ) with a high fraction of cold gas mass with respect to the stars. This fraction decreases very fast 
with decreasing redshift. Disk galaxies (blue) have a higher amount of cold gas at all times, especially at high redshift. 
Still, at low redshift there are many disk galaxies that have less than 10% cold gas. Those galaxies actually resemble 
SO properties and are not classical “disks.” 

In Table [8] we list the values of the cuts for the classification criteria. The values for the spheroids should be 
understood as an upper limit, while for the disks they represent a lower limit. 
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abs(circularity e) abs(circularity e) circularity e circularity e 

Fig. 19.— Left: cumulative fraction of all galaxies in dependency of the circularity e of the stars for all four redshifts. The colors reflect 
the different 6-value bins, as in Figure [5] From this we estimate the cut for the determination of the spheroidal part of a galaxy. This 
cut differs for each redshift. Right: cumulative fraction in order to determine the disk part of a galaxy. The values on the y-axis are the 
cumulative values from the corresponding point on the x-axis up to an e -value of 3. 



0 0.5 1 0 0.5 1 1.5 

^cold gas^star ^cold aas^star 



Fig. 20.— Left: fraction of the mass of the cold gas with respect to the stellar mass, both within the inner 10% of the virial radius. 
The cold gas fraction decreases continuously with decreasing redshift, for spheroids faster than for disk galaxies. Note that at z = 2 
most spheroidal galaxies have a mass fraction larger than 10%. At £ = 0.1 the disk galaxies are divided into two populations, one with 
a fraction higher than 30% and the other having less than 10%. We suggest that the second population comprises lenticular galaxies. 
Right: probability distribution function of the alignment of the galaxies compared to their DM halo in comparison with previous studies, 
as indicated in the plot. 


TABLE 8 

Upper/Lower Limits for the Classification Criteria 


Redshift 

£star 5 d 

£star 3 s 

Afcold/Afstar j d 

A/cold / -^star 3 S 

2 

0.55 

0.5 

0.5 

0.35 

1 

0.45 

0.55 

0.35 

0.2 

0.5 

0.45 

0.55 

0.275 

0.125 

0.1 

0.4 

0.6 

0.215 

0.065 


Note. — The values for the classification depending on redshift. We list the lower limit for the e s tar for disks ( d ) as well as spheroids 
(s). For the disks the M co \^/M s t ar is a lower limit, while the spheroids have to have less than the given value. 
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Fig. 21.— Left: A-parameter for the stellar component in the innermost 10% against that of the DM of the entire DM halo, at four 
different redshifts, divided into spheroidal (red circles) and disk (blue diamonds) galaxies. Right: A-distribution for the stellar component 
at four different redshifts, divided into spheroidal (red dashed) and disk (blue dot-dashed) galaxies. The stellar component of spheroids 
has a very low median value, whereas the disk galaxies have a higher one at all shown redshifts. 


TABLE 9 

K-S Test FOR THE A s tai-DlSTRIBUTIONS 


Redshift 

D 

Probability 

2 

0.558 

1.81 -10- 7 

1 

0.364 

3.42 -10 -5 

0.5 

0.443 

6.62 -10- 8 

0.1 

0.493 

2.40 TO -9 


Note. — The maximum distance D and the significance level (probability) resulting from the K-S test that the A s tar-distributions for 
disks and spheroids at different redshifts are drawn from the same distribution. 


APPENDIX B: THE ALIGNMENT OF THE GALAXY WITH ITS DM HALO 

The right panel of Figure [20] shows the probability distribution function of the cosine of the angles between the 
innermost 10% of the virial radius of our disk galaxies and the hosting DM halo. The solid blue line stands for the 
gaseous componen t, and the yellow dashed line for the stellar component. We extracted the other data points from 
ICodis ct al.1 ([20151) for an easier comparison with previous works. Overall, our results are in good agre e ment with 
previous stud i es. W e find that our disk galaxies h ave a slightly weak er alignment than th ose of lBett et al.1 (|2010H and 
iDeason et al.l (1201111 . Our results agree well with lHahn et all (12010(1 . ICodis et~ahl (12015H report slightly less aligned 
galaxies. 


APPENDIX C: THE SPIN PARAMETER A 

The left panels of Figure [21] show the A-parameter of the stellar component within the innermost 10% of the virial 
radius against that of the total DM halo for individual halos classified as spheroids (red circles) and disks (blue 
diamonds). We can clearly see the split-up of the two populations, where disk galaxies have a higher stellar spin than 
the spheroids. 

On the right side of Figure [5l] we show the A-distribution for the stellar component within the total virial radius. 
We also see a split-up of the spheroidal (red) and disk (blue) galaxies. The stars in disk galaxies have a higher spin 
than in spheroids. This seems plausible, since most stars in disk galaxies are found in the disk while most stars in 
spheroidal galaxies are found in the dominant bulge and thus are more spherically distributed. Hence, they have a net 
spin that is lower than that in disk galaxies, where the stars sum up their spin. This also shows that the stars trace 
the formation history of the galaxies. The stars in disks form out of the fresh gas with high angular momentum, while 
in spheroids many stars are accreted (via minor/major mergers). 

In Table [9] we show the results from the K-S test for the A s t a r-distribution. The two distributions for the disk and 
the spheroidal galaxies do not originate from the same distribution. 

In order to see how the A-parameter behaves for individual halos in the baryon and the DM-only runs, we identified 
the same halos in the baryonic and the DM-only run (as described earlier). The left panels of the left figure of Figure 
\77\ show the A-parameter of the DM component against the A-parameter of the DM-only run in logarithmic scaling 
for spheroids at redshifts z = 2 (upper panel) and z = 0.1 (lower panel). The values of the A-parameter of the 
corresponding halos in both runs are similar. The stellar mass within 10% of the virial radius of the baryon run is 
color-coded. On the right panels we find the AoM-values for the disks. At redshift z = 2 (upper panel) we note a slight 
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Fig. 22.— Left: A-parameter plotted logarithmically for the DM component, on the rr-axis the A for the baryon run and on the j/-axis 
for the DM-only run; the color codes the stellar mass within the inner 10% of the virial radius. On the left panels are the spheroids and 
on the right panels the disks. The upper panels show redshift z = 2 and the lower ones 2 = 0.1. Right: specific angular momentum (top 
panels) and the mass fractions (bottom panels) against distance from the center for the cold gas (blue), hot gas (red), stars (yellow), and 
DM (black) of our sample of galaxies at z = 0.1. On each left-hand side we show the analysis for the spheroidal galaxies, and on each 
right-hand side that for the disks. 



Fig. 23. — The radial profile of the 6-value averaged over spheroids (red circles) and disks (blue diamonds). 

tendency for the Adm to decrease wit h increasing stellar mass. The less massive disks tend to have higher A-values. 
This is in agreement with lBerta et all (120081) . wh o com p uted th e DM spin parameter of «52,000 disk galaxies from the 
SDSS (for details of the sample see references in iBerta et al.1 (120081) 1. They found a clear anticorrelation between the 
DM spin and the stellar mass, i.e., that galaxies with a low mass in general have higher DM spins. For the spheroidal 
galaxies we do not see this trend. In addition, this figure demonstrates that a high A does not automatically lead 
to a disk galaxy and a small A is no guarantee for a spheroidal galaxy. The scatter is equal for both galaxy types. 
The fact that Adm in the whole virial radius is similar in the DM-only and the baryon run is in good agreement with 
iBrvan et al.l (120131) . 


APPENDIX D: THE RADIAL SPECIFIC ANGULAR MOMENTUM PROFILES 

In order to illustrate the behavior of the specific angular momentum of the different components, we show on the 
right-hand side of Figure [25] the mean of our sample of spheroidal (left panel) and disk galaxies (right panel) at z = 0.1. 
In the upper panels we plot the differential specific angular momenta of the different components against radius, and 
in the lower panels we plot the individual mass fractions of the corresponding components against radius. The x-axis 
is linearly binned with a bin size of 10% of the virial radius of a halo. The specific angular momentum of the two 
galaxy types increases from the center of the halo up to the virial radius. On the outer parts the specific angular 
momentum of the hot gas follows the DM, while the specific angular momentum of the stars follows that of the cold 
gas. This might be due to the infalling substructures, which mostly consist of small gas clumps with stars that formed 
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out of it. The main difference of the two galaxy types becomes clear in the center: the angular momentum of the stars 
drops dramatically in spheroidal galaxies, while in disk galaxies it remains higher than that of the DM. With respect 
to the baryons, the outer region is dominated by the hot gas, which behaves similarly for spheroids and disks. 

Figure [23] shows the radial profile of the b- value. Therefore, we calculated 


b = log 10 


.7 star r )\ 

kpc km/s ) 



f Af star (< r)\ 

V ) ’ 


(Dl) 


averaged over the spheroids and disks. The slope is steeper for the spheroids than for the disk galaxies. This illustrates 
that the measurement of the kinematic properties is more sensitive to the radius for spheroidal than for disk galaxies. 
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